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1. Introduction 

The proton and the neutron are the fundamental building blocks of our world. They form 
heavier nuclei and thus the basis for atoms and are the only known source of stable baryonic mat- 
ter. Their structure can be studied by scattering leptons off nuclei and the most basic observables 
obtained from these processes are the nucleon form factors. The electromagnetic form factors are 
Lorentz-scalars which parametrize the matrix element of the electromagnetic current between two 
nucleon states at different momentum: 

{p'\qtq\p) = m)Fi{Q') + -^{{an)^aF2{Q^)^ (1.1) 

where ((^)) = u{p')^u{p) and = — = —{p' — p)^. rrif^ refers to the nucleon mass, and we 
only consider the isovector combination, i.e. the difference between u- and d-quark currents. In 
this case, contributions from disconnected diagrams cancel due to isospin symmetry. A different 
parametrization called the Sachs form factors, Ge{Q^) and Gm{Q^), is often used in the literature. 
We will discuss those in our upcoming paper |[l]]. 

In a similar way, the axial current can be parametrized in terms of two form factors. For the 
isovector form factors, i.e. the proton minus the neutron combination, they are called the axial form 
factor, Ga{Q^), and the pseudoscalar form factor, Gp{Q^): 

{p'\qysy^q\p) = {{fy,))GAQ^ + ^^^{{75))Gp{Q^). (1.2) 

The current work is based on mixed action calculations using two flavors of dynamical asqtad sea 
quarks and domain wall valence quarks. In previous years, we have reported on several other 
nucleon structure observables using this technology, see e.g. [gl, |, |]. We have studied form 
factors on full DWF lattices in Ref. [^] which also includes comparison to the work reported here. 
For a concise review of key results, we refer to Ref. [||]. For recent results from other groups, see 



Ref. n9n, and for recent reviews of the field Ref. [10, 11]. 



The current report focuses on selected results for nucleon form factors and uses several tech- 
nological updates. The final report using these improvements will be published soon [|l|] and will 
include several other major observables such as moments of generalized parton distributions and 
structure functions. 



2. Lattice technology 

As in our previous studies we employ the asqtad action for the sea quarks and the domain 
wall (DWF) action for the valence quarks. In addition, we also add one lighter mass to our data 
set. The tuning of the quarks masses and the choice of parameters have been discussed in Ref. [^]. 
The lattice spacing for all data sets corresponds to a = 0.124 fm with an uncertainty of 2%, see 
Ref. [|T|]. This yields a physical volume V = (2.5 fm)^ on the 20^ and V = (3.5 fm)^ on the 28^ 
lattices. The resulting physical values of the nucleon and pion masses are needed for our calculation 
and have been determined previously in Ref. Table [l] lists these numbers. 

In previous publications we often computed propagators by chopping each lattice in two halves 
and performing propagator calculations independently on both halves, cf. Ref. and references 
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Light mf^f 


Volume 




(am)N 


1 MeV 


mN / MeV 


0.007 


20^ X 64 


0.1842(7) 


0.696(7) 


292.99(111) 


1107.1(111) 


0.010 


28^ X 64 


0.2238(5) 


0.726(5) 


355.98(80) 


1154.8(80) 


0.010 


20^ X 64 


0.2238(5) 


0.726(5) 


355.98(80) 


1154.8(80) 


0.020 


20^ X 64 


0.3113(4) 


0.810(5) 


495.15(64) 


1288.4(80) 


0.030 


20^ X 64 


0.3752(5) 


0.878(5) 


596.79(80) 


1396.5(80) 



Table 1: Physical pion and nucleon masses. 



therein. In the present work we adopt a different technology and compute multiple source/sink pairs 
on a single gauge field. We find this approach both more convenient and more powerful, resulting 
in a superior statistical quality of our results. By choosing eight different source/sink pairs on a 
single gauge field, we managed to reduce our error bars by a factor of two. We also took possible 
sources of correlations into account by performing fits using the error correlation matrix among all 



data points on each ensemble, see e.g. |]6p, and the "super jackknife" technique, Refs. [ |13| , |14[ ], for 
combining data from different ensembles in a single fit. 



3. Form factor results 



We discuss several results of our calculation of nucleon form factors. To study the shape of 
the nucleon at large distances — which is a property that can be studied well by lattice calculations 
— we perform an expansion of the form factors at small Q^, yielding the mean square radii, (r?), 
as the slope, where / denotes either 1, 2, or A, corresponding to the form factor Fy, Fj, or Ga, 
respectively. Phenomenologically, the Dirac radius {r\) can be determined from fits to the form 
factor F\ {Q^) [ jlsl ] or from an analysis based on dispersion theory [ [T6| , [l7| , 18|. These two methods 
currently yield inconsistent results. For F2{Q^) there is a systematic discrepancy between spin- 
transfer and Rosenbluth experiments, the source is generally believed to be two-photon exchange 
processes, see Ref. [19] and references therein. Lattice calculations allow for a study without 
two-photon contamination and thus can be very useful in resolving this discrepancy. 

We have studied the form factors Fi{Q^) and 7^2 (2^) using dipole and tripole fits and also 
fit our lattice data to the simultaneous expansion in and mjt obtained from the small-scale 
expansion (SSE), see Refs. [20, 21] for the explicit form of these expressions. The advantage 
of the simultaneous expansion is that we do not make model-dependent assumptions on the 
dependence of the form factors, at variance with the use of dipole or tripole phenomenological 
formulae. The disadvantage is that the validity of the expansion will only hold for small values 
of and we have only few data points in that region. Thus, relying on the SSE expansion may 
increase the uncertainty, both statistical and systematic (which accounts for the unknown magnitude 
of higher order contributions). Since it is not feasible to determine all the low-energy constants 
involved in the chiral expressions by fitting to our lattice results, we fix some of them using their 
phenomenological values. 

We find that applying the cuts < 0.5 GeV^ and niji < 400 MeV yields an acceptable fit with 
X^/dof= 1.86 with two fit parameters. We could still describe the function well for larger Q^, but 
the dependence on m^t is worse. Despite this observation, we beheve that the apparent agreement 



3 



Nucleon form factors 



W. Schroers 



between the lattice data and the SSE form for > 0.5 GeV^ is merely accidental since we have no 
reason to believe that the SSE at the order given is valid for that range of Q^. Figure |l] summarizes 
our results. The left panel shows the isovector form factor F\ {Q^) lattice data with the best fit SSE 
at triji = 292.99 MeV. The right panel shows the resulting chiral extrapolation of the Dirac radii as 
a function of the pion mass, mji. For illustration purposes, we have also included the Dirac radii 
obtained from dipole fits in the graph. However, these data points have no influence on the curve 
presented and simply compare the two fit methodologies. The red star shows the empirical value 



taken from Ref. |15]. 




0.5 0.75 
d /Ge\/ 



1.5 



0.5 




Figure 1: Isovector form factor Fi (g^) lattice data with best fit SSE at nin = 292.99 MeV (left panel). The 
line in the right-hand panel shows the resulting Dirac radii, (r^ ). Also shown as the data points are the Dirac 
radii obtained from dipole fits to the form factors at each pion mass 



The corresponding fit for the isovector F2{Q^) had a quality of ;^^/dof= 1.31 with four fit 
parameters when applying the cuts < 0.3 GeV^ and niji < 400 MeV. However, we noticed that 
the dependence of the SSE expression was in worse agreement than in the case of Fi{Q^), 
while the niji dependence was better. The resulting Pauli radii and anomalous magnetic moments 
are shown in Fig. The left panel shows the Pauli radius (rj), while the right panel shows the 
isovector anomalous magnetic moment, Ky. Again, data points from tripole fits are included in 
the plot, but have no influence on the fit. We find that (r^) is described well by the fit and the fit 




100 200 300 400 500 600 700 100 200 300 400 500 600 700 

m^/MeV m^/MeV 



Figure 2: Isovector Pauli radius {r^} (left panel) and isovector anomalous magnetic moment ff,, (right panel) 
as a function of the pion mass. 



even gets close to the experimental point. The magnetic moment Ky is lower than the experimental 
value. We find these fits encouraging, but believe that the excellent agreement with experiment is 
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accidental for {r\). Data with smaller pion mass are needed to verify this finding since K^. does not 
yet agree, although it originates from the same data set. 

Similar to our fit strategy for the vector form factors, we also adopt a simultaneous fit to the 

and /M;r-dependence of the axial form factor, Ga{Q^)- Figure ^ shows the result of the chiral fit 
together with the dipole fit and the data set for the 28^ lattice at niji = 355.98 MeV with a fitting 
range of < 0.4 GeV^ for the SSE expansion and all values for the dipole fit. The SSE fit 
gives a ;^^/dof= 1.73. The resulting axial radius, however, is (r^) = 0.1560(60) fm^ for the cuts 

< 0.4 GeV^ and mji < 400 MeV, substantially lower than the experimental value in Ref. [^]. 
Future lattice calculations at smaller pion masses will be crucial to resolve this issue. 



— Dipole fit 

- SSE tit 



"""0 0.1 0.2 0.3 0.4 0.5 0.6 
cf/GeV^ 

Figure 3: Comparison of dipole and chiral fit to Ga{Q^) for the 28^ lattice at m^i = 355.98 MeV. 

The induced pseudoscalar form factor, Gp{Q^), is not described by a dipole type fit formula. 
Instead, it is commonly fit using a pion-pole expression, giving excellent agreement with the data, 
see Ref. [p^]. We performed two kinds of fits: First, we repeated the analysis done previously with 
the other form factors, i.e. performing a combined fit in and mj^. Second, we took the pion-pole 
form as a function of and fit it using a single ensemble with fixed m-n, treating the pion mass as 
a free parameter. 

In the first case, we again find that kinematic cuts of < 0.5 GeV^ and niji < 400 MeV yield 
reasonable results. In the second case, we find the location of the pion pole, nijc = 417(43) MeV, 
with an uncertainty of 10% within the actual value of niji = 355.98 MeV on the 28^ lattice. The 
resulting ;^^/dof= 0.94 indicates an excellent fit to the data. Figure ^ shows a comparison of the 
two fits. It is evident that both fits manage to describe the data well, but the uncertainty of the curve 
with as a free parameter is larger. We conclude that the isovector Gp{Q^) is described well by 
the pion-pole form. 

4. Summary and Outlook 

We have computed the isovector form factors of the vector and the axial current for the nucleon 
within our framework of mixed action calculations. We have applied new analysis techniques that 
result in substantially reduced error bars at minimal additional computational cost. Furthermore, 
due to the agreement between the lattice results reported in this proceeding and the corresponding 
one reported in |^ obtained using full DWF, we are confident that the hybrid calculations employed 
do not suffer from systematic effects. 
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Figure 4: Comparison of pion-pole fits to isovector Gp{Q^) on the 28^ lattice with fixed pion pole and with 
the pion pole as a free parameter. 



We find that a combination of chiral fits and lattice data is possible with the cuirent generation 
of lattice calculations. This way, we obtain qualitative agreement with many features we expect to 
hold when approaching the chiral limit. We expect that the upcoming generation of lattice calcu- 
lations will provide conclusive quantitative results from first principles that will be in agreement 
with experiment without resorting to assumptions on functional behavior outside of what can be 
predicted by chiral perturbation theory. We are able to provide fits to the vector form factors, 
Pi{Q^) and F2{Q^). While we have no explanation for the discrepancy of the axial radius, {r\), we 
find that induced pseudoscalar form factor is described well by the pion-pole form. 
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